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^ ' Abstract: Acoustic peaks in the spectrum of the cosmic microwave background in spher- 

^1 ■ ically symmetric inhomogeneous cosmological models are studied. At the photon-baryon 

^^ ■ decoupHng epoch, the universe may be assumed to be dominated by non-relativistic mat- 

o 



ter, and thus we may treat radiation as a test field in the universe filled with dust which 
is described by the Lemaitre-Tolman-Bondi (LTB) solution. First, we give an LTB model 

f^ I whose distance-redshift relation agrees with that of the concordance ACDM model in the 

whole redshift domain and which is well approximated by the Einstein-de Sitter universe 
at and before decoupling. We determine the decoupling epoch in this LTB universe by 

^ ' Gamow's criterion and then calculate the positions of acoustic peaks. Thus obtained re- 

j^ ■ suits are not consistent with the WMAP data. However, we find that one can fit the peak 

positions by appropriately modifying the LTB model, namely, by allowing the deviation of 
the distance-redshift relation from that of the concordance ACDM model at z > 2 where 
no observational data are available at present. Thus there is still a possibility of explaining 
the apparent accelerated expansion of the universe by inhomogeneity without resorting to 
dark energy if we abandon the Copernican principle. Even if we do not take this extreme 
attitude, it also suggests that local, isotropic inhomogeneities around us may seriously af- 
fect the determination of the density contents of the universe unless the possible existence 
of such inhomogeneities is properly taken into account. 
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1. Introduction 

In the study of cosmology, we usually assume that we do not live in a special position in the 
universe, the so-called Copernican principle. However, although this principle is natural, 
it should not be blindly assumed but should be justified observationally. Conventionally 
observational data are interpreted under the assumption of a homogeneous and isotropic 
universe on average. Therefore it is not clear at all how big the systematic errors would be 
in the determination of the cosmological parameters if this assumption were abandoned. 
In other words, it is important to investigate possible "anti-Copernican" models of the 
universe and test if such models can be observationally excluded. In this paper, as one 
of such attempts, we try to construct a cosmological model without dark energy which is 
consistent with the observed distance-redshift relation as well as with the WMAP data. 

In anti-Copernican models, we are assumed to be located at a special place in the 
universe, usually at the center of a spherically symmetric inhomogeneous universe [0, ||, |3|, 
Q, |5|. In recent years, such models have attracted much attention |^, |^ ^, ^, ^, |ll], |l2|, 

0, g, 0, 0, 0, |i8|, 0, ig, 0, ig, ^ 

,|2 



27, 25 



31 



4 35, 36 



, 24|, and various ways to observationally test these 
models have been proposed by many authors 

1, 0,11, 11,11, 11,11,13, HI. 



33, ^, |39| 

One of the simplest ways to construct an anti-Copernican model is to solve the inverse 
problem of the distance-redshift relation. Although the isotropy of the universe around 



us has been confirmed with high accuracy by the observation of the cosmic microwave 
background (CMB), this does not automaticahy imply homogeneity of the universe. Thus, 
in solving the inverse problem, we may assume that the universe is spherically symmetric 
around us. In addition, we usually assume that the universe is dominated by cold dark 
matter, that is, by dust. The spherically symmetric dust filled spacetime is described by the 
Lemaitre-Tolman-Bondi (LTB) solution. The LTB solution has three arbitrary functions 
of the radial coordinate approximately corresponding to the density profile, the spatial 
curvature and the big-bang time perturbation, with one of them being a gauge degree of 
freedom representing the choice of the radial coordinate. These arbitrary functions may 
be determined by requiring that the resulting LTB universe be consistent with selected 
important observational data (e.g., the distance-redshift relation). However, we should 
note that it is not apparent at all if these three functions have enough degrees of freedom 
to fit all of the important observational data. 

In 1999, Celerier solved the inverse problem analytically at small redshifts z <C 1 
in the form of the Maclaurin series [^. Then, in 2002, Iguchi, Nakamura and Nakao 
constructed numerically an LTB model whose distance redshift relation agrees with that 
of the ACDM model at z < 1.6 [^] ^. However, they could not go beyond z ~ 1.6 due 
to a technical problem. Since then the inverse problem has been discussed by various 



authors |51 



55, 5C, 57 1 . In 2008, Yoo, Kai and Nakao succeeded in constructing 
an LTB model whose distance-redshift relation agrees with that of the concordance ACDM 
model in the whole redshift domain and which is homogeneous in the early stage of the 
universe, that is, the big-bang time perturbation being taken to be zero p^. Note that 
because of the existence of a gauge degree of freedom in the choice of the radial coordinate. 



there remains only one functional degree of freedom. In [57|, this is represented by the 
function that determines the spatial curvature. 

Not only the distance-redshift relation but also CMB observations in inhomogeneous 
cosmology have been widely discussed 0, |, |, 0, |8|, |3|, 0, |l|, |l|, ||, g^, |5|, 0, |2|, ||, ||] . 
In many of these works, it is assumed that the universe is homogeneous in the spacelike 
asymptotic region from which the CMB photons come, hence the CMB photon distribution 
is assumed to be the same as that in the homogeneous and isotropic universe models at 
the time of decoupling. Several authors proposed parametrized LTB models and gave 
constraints on the parameters with observational data P, ^, 10, 11, 14, |^, 16, ^]. 



Recently, Bolejko and Wyithe |l^ suggested that it is possible to construct an LTB 
model which is consistent not only with the distance-redshift relation but also with the 
acoustic peak positions of the WMAP data [pO| by choosing appropriately the arbitrary 



functions of the LTB solution. In this paper, by modifying the LTB model given in Ref. |57| 
we show explicitly that this is indeed possible. For simplicity, we use Gamow's criterion to 
determine the decoupling epoch instead of invoking a precise numerical calculation. This 
is the same procedure adopted in Ref. [^2| . An advantage of this simplified prescription is 
that it makes it easy to understand the physical degrees of freedom in the LTB universe 
that determine the CMB anisotropy spectrum. We note that, in our construction, only the 



^Before Refs. 



49[, Mustapha et.al discussed an inverse problem using the LTB solution fcol 



asymptotic homogeneity on the past hght cone of the observer at the symmetry center was 
necessary in contrast to a much stronger asymptotic spatial homogeneity condition (see a 
related discussion in Ref. ||6l[|). 

This paper is organized as follows. In §Q, we give a brief review of the LTB solution 
and give a fitting function for the spatial curvature of the LTB model that has the exactly 
same distance-redshift relation as the concordance ACDM model [^. In §^ we discuss 
Gamow's criterion and the physical degrees of freedom at decoupling. The position of the 
first acoustic peak in the CMB spectrum in an LTB universe is discussed in §^. In §|5|, 
we construct an LTB model which is consistent with the acoustic peak positions of the 
WMAP data by modifying the asymptotic structure of the LTB model obtained in |p7|| . 
§^ is devoted to summary and discussion. 

2. LTB model from the inverse problem 

As mentioned in the introduction, we consider a spherically symmetric inhomogeneous 
universe filled with dust. This universe is described by an exact solution of the Einstein 
equations, known as the Lemaitre-Tolman-Bondi (LTB) solution. The metric of the LTB 
solution is given by 

1 — k[r)r'^ 

where k{r) is an arbitrary function of the radial coordinate r. The matter is dust whose 
stress-energy tensor is given by 

T''" = pu^u\ (2.2) 

where p = p{t,r) is the mass density, and u"" is the four-velocity of the fiuid element. The 
coordinate system in Eq. ( |2.1| ) is chosen in such a way that u^ = (1,0, 0, 0). 
The area radius R{t, r) satisfies one of the Einstein equations, 

(2Ey=^jm^t(ry, (2.3) 



\dt J R 

where M(r) is an arbitrary function related to the mass density p by 



/ X 1 dM(r) , ,, 



Following Ref. |62], we write the solution of Eq. ( |2.3| ) in the form, 

R(t, r) = (6GM(r))i/3(t _ t^(r)fl^S{x), (2.5) 

x = ^k(ry(^-^^^f^\ (2.6) 



V 6GM(r) 

where t^ij^ is an arbitrary function which determines the big bang time, and S{x) is a 
function defined implicitly as 



Six) 



cosh \/ —r\ — 1 — (sinh -v/^r? — \/— r?)^'^ ^ 

\ x = -^ — ^7-—^^ for x < . 



6-*^/^ (sinh -^—r/ — ^— r/)^/^ 6^/^ 

1 — cos y^ (y^ — sin ^Jfl)^l'^ 

, 6V3(^_sin^)2/3' ^= ^^73 



(2.7) 
for X > , 



and jS'(O) = (3/4)"'^'^. The function S{x) is analytic for x < (vr/S)^'^. Some characteristics 



of the function S{x) are given in Refs. [^] and [62]. 

As shown in the above, the LTB solution has three arbitrary functions, k{r), M{r) 
and tB{r). One of them is a gauge degree of freedom for rescaling of the radial coordinate 
r. In this paper, we fix this by setting 

M{r) = ^npor^ (2.8) 

where po is the energy density at the symmetry center at present po = p{tQ,0). As in 
the case of the homogeneous and isotropic universe, the present Hubble parameter Hq is 
related to po as 

H^ + k{0)c^ = ^vrGpo. (2.9) 



As in Ref. |57|, we assume the simultaneous big bang, i.e., 

tB{r) = 0. (2.10) 

For notational simplicity, we introduce dimensionless quantities, 

Hor ~ k{r)c^ 

r := , k{r) : — 



/72 



The observed distance-redshift relation is consistent with the homogeneous and isotropic 
universe model with (r^mOj^Ao) = (0.3,0.7), where Qmo is the density parameter of the 
total non-relativistic matter (i.e., cold dark matter plus baryons) and Oao is that of the 
cosmological constant, the so-called concordance ACDM model. We determine k{f) so that 
the distance-redshift relation of our LTB model agrees with that of the concordance ACDM 
model. In Ref. ||57|] , the inverse problem was solved numerically for k{f). A fitting function 



to the numerical result obtained in Ref. |57] is given by'^ 

0.545745 2.22881 

0.211472 + Vo:026176 + f ~ (0.807782 + ^0^026176 + f ) 



^fit(^) = ^^,.,^^ ; /77^^^,^. , ,, - • ^ . . (2.11) 



As shown in Fig.|l], the distance in our LTB universe model with k{f) = k^x.i'f') agrees with 
that in the concordance ACDM model in the whole redshift domain. 

3. Decoupling Epoch in the LTB universe 

The decoupling between photons and baryons occurs in an inhomogeneous universe just 
as in the case of a homogeneous universe. That is, it occurs when the mean free path of 



^It should be noted that the fitting function has a non- vanishing first derivative at the center, corre- 
sponding to the existence of a spike in the density profile. In our model the central region has a void-like 
structure with f^mo ~ 0.09 at the center. However, the density quickly rises to SlmO ~ 0.2 a.t z = 0.3 in 
agreement with large scale structure observations, as shown in Ref. [p7|]. The longitudinal Hubble parameter 
Hl := dtdrR/drR in our LTB model is about 10% larger than that in the concordance ACDM model at 
z ~ 2. Although we need more careful investigations, considering large error bars and uncertainties in the 
local Hubble measurements |63|, there seems no apparent confiict between our model and observation. 
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Figure 1: Distance modulus for our LTB model and the concordance ACDM model. 



photons becomes effectively infinite due to almost complete recombination of electrons to 
protons. 

Since there is no radiation component in our LTB model, we cannot treat the decou- 
pling in a rigorous manner. However, as in the concordance model, we expect the energy 
density of the radiation to be only a small fraction of the total density at decoupling, hence 
its effect on the spacetime geometry is small, if not negligible. In fact, the radiation energy 
density estimated in our LTB model turns out to be about 20% of the total energy density. 
This means that treating the radiation as a test field in our model is consistent to a first 
approximation. 

Another approximation we adopt is the instantaneous decoupling Namely, we assume 
decoupling to occur on a single spacelike hypersurface. Since our LTB universe model 
is inhomogeneous but spherically symmetric, it is natural to assume that the decoupling 
hypersurface is also inohomogeneous but spherically symmetric. Thus it is specified by the 
form, 



t = iD(r) . 



(3.1) 



The cross section of this hypersurface with the past directed null cone from the observer at 
the center constitutes the last scattering surface (LSS) of CMB photons (see Fig. |2|). The 
LSS is a spacelike 2-dimensional sphere by the assumed symmetry. We use the subscript * 
to express physical quantities on the LSS. In our approximation, ignoring secondary effects, 
the CMB anisotropy is essentially determined by the distribution of photons on the LSS. 

The geodesic equations to determine the past light cone from the observer at the center 



time 




Decoupling epoch 
(hyper surface) 



LSS 
(sphere) 



Figure 2: Schematic figure for tire liypersurface of tlie decoupling epocli in an LTB universe. 



are written in the form, 



[l + z] 



dt 



dr-R 



dz dtdfR 



dr C\l\ — k{r)r'^ 



(3.2) 
(3.3) 



where the past directed radial null geodesies have been parametrized by the cosmological 
redshift z. We denote the solution of the above equations by 



t = tic{z) , r = ric{z). 



(3.4) 



Now to discuss the decoupling condition, we consider the matter contents of the uni- 
verse around the time of decoupling. For simplicity, we consider a universe consists of cold 
dark matter, protons and electrons, and neutral hydrogen atoms. In particular, we neglect 
helium. Since the contribution of these other components is not large, this simplification 
should not lead to a serious error in our analysis. We also assume that, until the decoupling 
time, photons, electrons and protons are in thermal equilibrium. The energy density of 
the electrons and photons is negligible, and hence the constituents of our LTB model are 
cold dark dark matter and baryons, the latter of which consist of protons and hydrogen 
atoms. Thus the baryon number density n-b is equal to the total number density of protons 
and hydrogen atoms, and the electron number density Ue is equal to the proton number 
density. 



In homogeneous and isotropic cosmology, the decouphng time is well determined by 
Gamow's criterion, 

H = r, (3.5) 

where H is the Hubble parameter and F is the rate of collisions of a photon with electrons. 
Using the Thomson scattering cross section ctt and the electron number density Ug, T is 
written as 

r = creeO"T • (3-6) 

In our LTB model, we also adopt this Gamow's criterion, with the identification of the "Hubble 
parameter H" with 

i/2 = __p. (3.7) 

Since, by virtue of Eq. (2.10), our LTB universe model will be very similar to the Einstein- 



de Sitter universe near the decoupling time, the above definition should be accurate enough 
for the present purpose. 

In order to estimate the electron number density, let us consider the ionization rate 
Xf. := neju^,- In the thermal equilibrium, the ionization rate Xe satisfies Saha's equation, 

1-Xc 4^/2C(3) fkBTV^^ /l3^9eV\ 

where C{x) is the zeta function, and T, r/, /cb and me are the temperature, the baryon- 
to-photon ratio, the Boltzmann constant and the electron mass, respectively. Since the 
ionization rate at decoupling drops down to Xg ~ 10~^, we approximate the above equation 

by 

Using this equation, Gamow's criterion ( |3.5[ ) is rewritten in the form, 

32^/2^C(3) Gp /^bTo y/' /To^/' /13^59eV\ ,^ ^^, 

where n^o and Tq ~ 2.725K are the present photon number density and observed CMB 
temperature, respectively. 

Since we assume thermal equilibrium of electrons, protons and photons until the decou- 
pling time t = t]:,{r), the physical state of the decoupling hypersurface, which is spherically 
symmetric, is determined by the distributions of the temperature T = To^r), the baryon- 
to-photon ratio rj = r/D(r), and the matter energy density p = p{ti:){r),r). For convenience, 
in place of p{tD(r),r), we introduce the following quantity: 

p(tD(r),r) f To y 
The quantity aB{r) is proportional to the ratio of the matter density and the photon 



number density. Then, from Eq. (3.10), we obtain 

3/2 / T" \ 3/2 



32V2^C(3) Gpo fkBny f^f'^oY' f U.59eY \ 



We note that once 7b (r) and ai:){r) are given, the hypersurface t = tD{r) can be obtained 
from Eq. (|3lll ). 



On the LSS, we must have 



1 + z*. (3.13) 



To To 

Then, if we regard Tb(r), r]Y)(r) and aol^) as mutually independent functions, we have 



one functional condition ( |3.12| ) and one boundary condition at z = z^:. to constrain these 
three functions. But these are not enough to determine the decoupling hypersurface, t = 
t]:,{r), through Eq. (3.11). We would need two more functional conditions. However, 



for the present purpose, we do not need to determine the total shape of the decoupling 
hypersurface, but we only have to know the location of the LSS. Therefore we only need 



a single condition on the LSS in addition to Eqs. (|3.12D and ([3.131 ). In this paper, for 
simplicity, we impose 

77* = 6.2 X 10"^° . (3.14) 

We numerically solve the radial null geodesic equations ( |3.2D and ( |3.3[ ) . At each redshift 
z, we define T and a by 

r = (i + z)ro, (3.15) 

During the numerical integration, we check at each redshift z whether Eq. ( [3. 121) is satisfied 



by Td = r, aD = ce and 77D = ?7*. If Eq. (3.12) is satisfied, we stop integrating Eqs. (p^) 



and (3^), and identify z, T and a at this moment with z*, T* and a*, respectively. The 
result we have obtained is 

-f = l + z*~1129, (3.17) 

Jo 

a=,~ 4.856, (3.18) 

for our LTB universe model. 

In the above discussion, we have focused on only a* and r/*, and we have given no 
further constraint on a-Q{r) and r]D{r). Therefore, we can assume a]:){r) = a* and 7?d(^) = 
r]^: over the whole decoupling hypersurface. Inhomogeneities of qd and ijd correspond to an 
isocurvature perturbation of the matter density and a baryon isocurvature perturbation, 
respectively. However this assumption might lead to a contradiction to the observational 
data of the kinematic Sunyaev-Zel'dovich effects as will be discussed in a forthcoming paper. 
It is also worthy to notice that the anomalous primordial Lithium abundances reported in 



Ref.|64] can be explained by introducing the baryon isocurvature perturbation p^ . 



4. acoustic peaks in the CMB spectrum 

4.1 acoustic peaks in the homogeneous and isotropic background 

Let us first briefiy review the positions of the acoustic peaks in the CMB spectrum in the 
homogeneous and isotropic background universe. The acoustic peak positions in the CMB 



spectrum can be written as fe 



Im = {m - (t)m)T^ , (4.1) 



' s* 



where d^ is the angular diameter distance from the LSS to the observer at present, rg* is 
the radius of the sound horizon at the time of decoupling, and ipm is a small correction to 
the position of the TTi-th peak. As explicitly shown in the Appendix ^, r^^ and (f)m (for 
m = 1, 2, 3) can be expressed in terms of Wm := ^mo^^j '^b '■= ^bo^^j ^-t '■= ^-yoh?, z^^ and 
the spectral index Ug which characterizes the spectrum of the initial density perturbation, 
where ilmOi ^bO; ^^70 ^-re the density parameters of the total non-relativistic matter, of 
baryons and of radiation, respectively, and h = ifo/(100km/s/Mpc). 

In order to compare the first acoustic peak in the CMB spectrum predicted by our 
LTB universe model to the observational results, it is convenient to use the quantity S 
defined by [^] 

"^ •= ^WMAP ' (^•^) 

where ii is the first acoustic peak in the CMB spectrum and ^wmap _ 220.8 is the mean 
value of the first peak in the WMAP data ||6^ (see Table |l]). Since the la error in the 
WMAP data is Aii = 0.7, models with |5 - 1| > 0.0063 are ruled out at 67% CL. 

Before closing this subsection, it may be worthwhile to note the role of the instan- 
taneous decoupling approximation. Its role is to identify a homogeneous universe model 
which is a good approximation to our universe model around the time of decoupling. Once 
the parameters of a homogeneous universe model are identified, the positions of the acous- 
tic peaks are calculated from the fitting formula (30) derived from numerical simulations 
for homogeneous universe models. Hence except for the effect of errors in the estimation 
of the cosmological parameters, there will be no deviation of the acoustic peaks due to the 
instantaneous decoupling approximation. 

4.2 Position of the first acoustic peak in the LTB model 

Because the inhomogeneity in our LTB model consists of only growing modes, the inho- 
mogeneity decreases as we go back to the early stage. Actually, our LTB model can be 
well approximated by the Einstein-de Sitter (EdS) universe at the decoupling time. In 
this subsection, we therefore approximate our LTB model at the decoupling epoch by an 
EdS universe model. Hereafter, we place a bar over a quantity in this EdS universe. We 
determine it by setting the decoupling temperature To and the baryon-to-photon ratio f] 
in this EdS universe model equal to T* on the LSS in our LTB model and 77* given by 
Eq. ( |3.14 ), respectively. Further, we set the present CMB temperature to be equal to To 
also in the EdS model. These three conditions determine O^, oJb and Cj^, and the Hubble 
parameter ^0 at present, 

j-2 _ zt2 f'^o\ _ 8vrG [To 



^»-'^.is^"f'- s) ■ ("' 



Using the definition of a, we have 



Hi = ^poa* = [1 + fe(0)]i/o^a, ~ OmUOH^a^ . (4.4) 



Since the present CMB temperature is Tq, we have 



^0 A 



\ 4 



tj^ = 4.2 X 10"^ ^ . (4.5) 

^ V2-725Ky ^ ^ 

The quantity tDb is expressed in terms of the baryon-to-photon ratio as 

8ttG _ jo 7 f To \^ , , 

Wb = ■^^r]*'n'yomph = 3.658 x 10 r/* I 2 725K j ' ^ ^^^ 

where ?Tip is the proton mass. As for the matter density, we have Q^q = 1 by assumption, 
hence 

oj^ = n^oh'^ = h^ = OmUOa^h'^ , (4.7) 

where we have used Eq. ([4.41) . 

By using Eqs. ( |3.14| ) and (|3.18| ), we can evaluate uij, Wb and (Dm- Then, substituting 
Lo-y = ui^, Wb = (Db and LOm = Wm for Eqs. ( A.ll| ), ( |A.12| ), ( [A. 15 ) and ( A.16| ), we find 



|5 - 1| ~ 0.075 > 0.0063, (4.8) 



where we have used Eq. (3.17) and set h = 0.71 and To = 2.725K. This result imphes that 



our LTB model cannot explain the observed first peak position in the CMB spectrum. 

5. Matching peaks in the CMB spectrum 

Let us consider if we can resolve the inconsistency of our LTB model with the observed 
first acoustic peak position in the CMB spectrum. Since SNe observations of the distance- 
redshift relation are still restricted to relatively low redshifts z < 2, the distance-redshift 
relation may not necessarily agree with that in the concordance ACDM model at z > 2. 
Therefore we may modify the curvature function k{f) in the domain of large r so that the 
position of the first peak in an LTB model agrees with the observed one without affecting 
the distance-redshift relation at 2: < 2. 

As soon as we allow modifications of our LTB model at high redshifts, we have too 
much freedom to fix the model uniquely. Therefore, for simplicity, we consider a single- 
parameter modification. Namely, we adopt the curvature function, 

~k{f;A) = ht{r)xf{r;A), (5.1) 

where the function /(f; A) is defined by 



fir; A) 



1 for f < 2 , 

16A(r-2f (323-123r + 12r'^) , .^ „^ 

1 + ^ ^ \,,, ^ for 2<f<9/2, (5.2) 

3125 
1 + A for 9/2 < f . 



10 
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Figure 3: fc(f; A) and f{f; A) as functions of f for several values of A. 



The functions k{f; A) and /(f; ^) are depicted in Fig. |3| for several values of A. 

As in the previous section, we adopt the value rj^ = 6.2 x 10^^*^. Then we can evaluate 
5 — 1 for each value of A. 5 — 1 is depicted as a function of A in Fig.H. As shown in 
this figure, S almost vanishes at yl = —1.069. This fact means that the modified LTB 
model with t^ = 0, k = k{f; —1.069), 77* = 6.2 x 10^^'^ and h = 0.71 is consistent with the 
observed distance-redshift relation and the position of the first acoustic peak in the CMB 
spectrum simultaneously. Then we may check if the positions of the second and third peaks 
are also consistent. As shown in the Table |l|, they indeed turn out to be consistent with 



11 
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Figure 4: S — l(left), z*(center) and ft,(right) as a function of A. 



the WMAP data |6^. The distance modulus in this LTB universe model with A = —1.069 
is shown in Fig. |5|. 

The resultant LTB universe model is approximated by an EdS universe model with 
the Hubble parameter h ~ 0.49 in the spatially asymptotic region. This result is consistent 
with previous work in the literature [pq , |6|, |9|, |l^ in which the possibility for lower Hubble 
universes to explain the CMB spectrum without dark energy has been pointed out. 
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Table 1: second and third peak positions in WMAP data and our LTB model with A 
We have assumed «= — 0.963. 
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Figure 5: The distance modulus for the LTB universe model with t^ = 0, fc(f) 
^fit(r)/(r; -1.069). 



6. Summary and Discussion 

In this paper, we investigated if an LTB universe model can account not only for the 
observed distance-redshift relation but also for the observed peak positions in the CMB 
spectrum, assuming that we are at the symmetry center of the universe. 

First, we presented an LTB model whose distance-redshift relation is equal to that of 
the concordance ACDM model in the whole redshift domain. This model is expressed in 
terms of a fitting function for the curvature numerically obtained in Ref. [57|. Then, we 
determined the last scattering surface for photons in this model by numerically integrating 
the past-directed radial null geodesies emanating from the central observer until Gamow's 
criterion is satisfied. 

The metric of this LTB model can be well approximated by the EdS universe at the 
decoupling epoch. Thus, assuming that our model universe experiences the same history as 
that of a homogeneous and isotropic universe before decoupling, we evaluated the positions 
of the acoustic peaks in the CMB spectrum. It was found that thus obtained position of the 
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first acoustic peak deviates from the observed position more than 7%, implying that our 
LTB universe model whose distance-redshift relation agrees with that of the concordance 
ACDM model for the entire redshift domain is ruled out. 

In order to resolve this inconsistency, we considered modifications of the LTB model 
in the region far from the symmetry center, corresponding to the redshift domain of z > 2. 
Since the observation of Type la supernovae are limited within z < 2, this modification 
does not contradict the current observational data. Then we considered a one-parameter 
family of LTB models whose distance-redshift relation matches the observational data at 
z < 2, and found a parameter that gives the first acoustic peak position consistent with the 
observed peak position. Then we checked the second and third peak positions and found 
that they are also consistent with the WMAP data |60|1 . 

We have not discussed relative heights of acoustic peaks in this paper. At the moment, 
we have no idea if one can also fit them by an LTB model. There is, however, a hint in the 



literature. In Ref. |66|, it was reported that the observed CMB spectrum can be regarded 
as that in the EdS universe with a low value of the Hubble constant if one tunes the 
primordial power spectrum of the density perturbation. If this is the case, then we may 
be able to find an LTB model which can account for the observed relative heights of the 
acoustic peaks. This issue is left for future work. 

In this paper, we have not introduced isocurvature components in the inhomogeneity 
on the last scattering surface. In a forthcoming paper, we plan to investigate if we may con- 
struct an LTB model consistent with the observation of the kinematic Sunyaev-Zel'dovich 
effect. Our preliminary analysis indicates that isocurvature perturbations between the 
non-relativistic matter and photons are necessary in order to explain it. However, if we 
introduce large radial isocurvature perturbations on the decoupling hypersurface, the his- 
tory of the inhomogeneous universe model before decoupling would differ substantially 
from the conventional adiabatic perturbation scenario in the homogeneous and isotropic 
background. In such a case, our analysis of the acoustic peaks in the CMB spectrum may 
be invalidated. However, since we have considered only a single-parameter family of LTB 
models, while there still remains a functional degree of freedom, it is premature to make 
any definite statement at the moment. Apparently much more work seems necessary to 
exclude or select the LTB universe as a viable alternative model of our universe. 
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A. peak positions in homogeneous universes 

A.l expression for the sound horizon 

Since the universe considered in this paper is well approximated by the homogeneous and 
isotropic universe with vanishing spatial curvature near the decoupling epoch, we consider 
the only such a universe model here. The infinitesimal world interval is given by 

ds^ = -dt^ + a^{t){dr^ + r^dn^). (A.l) 

The scale factor a{t) is normalized so that its present value gq is unity. 
The sound horizon rg is defined by 

/'* -^' = r -^u 

/o a Jq a'^H 

where Cg is the sound speed given by 

c 



rs= / -dt= I -:^da, (A.2) 



VH^ + rY 



(A.3) 



Ap^ Auj^ T Auj^ I + z 
At the matter-radiation equality time, we have 

^e, := p^^ = 1^, (A.5) 

where we have used that the temperature T^q at the equality time can be written as 

Teq = ^To. (A.6) 

UJ-y 

The denominator of the integrand in ( A.2| ) can be written as 



a'^H = Jnn,oacq(^ + l], (A.7) 

where Oeq is the scale factor at the equality time, Omo is the density parameter of the total 
non-relativistic matter and we have used the relation R = aReq/aeq- Using the expressions 
(|A.3|) and (|A.7|), we can perform the integration in ([A.2|) as 



2 / V 1 + i?* + y -R* + -Req 



In ' : %^ "" , (A.8) 



fceq V 3/Leq \ 1 + \/ Ri 



eq y "-'-' "'eq \ i. -r y i tcq 



where 



and 



= ^ /2fVo ^ ^ U^^o^ = 4.714 X lO-^^/^V^Mpc-i , (A.9) 

R^ = pL^, (A.IO) 
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A. 2 small shift parameters for the first, second and third peaks 

We define the ratio r^ of the radiation density to the matter density at decouplin, 



Pn 






UJr, 



-(1+^*). 



Using this ratio, a fitting formula given in Ref. |67] gives 



where 



(1)1 + 5> 



eirl^ 



0.037(ns - 1), 



C2r-^3+o.05(ns-l), 



ai = 0.286 + 0.626a;b, 

02 = 0.1786 - 6.308a;b + 174.9a;g - 1168wb, 
Co = -0.1 + 0.213exp(-52wb), 
0.063 exp(-3500w^) + 0.015, 



ci 



C2 = 6 X 10^'^ + 0.137(wb - 0.07)^ 

C3 = 0.8 + 70wb, 

ei = 0.302 - 2.112wb + 0.15exp(-384wb), 

62 = -0.04 - 4.5a;b. 



g as 



(A.ll) 



(A.12) 
(A.13) 

(A.14) 



(A.15) 
(A.16) 
(A.17) 
(A.18) 
(A.19) 
(A.20) 
(A.21) 
(A.22) 
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